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Abstract. I describe the advantage of using singular value decomposition as a 
diagnostic tool for exploring the potential and limitations of seismic data. Using 
stellar models coupled with the expected errors in seismic and complementary 
data we can predict the precision in the stellar parameters. This in turn allows us 
to quantify if and to what extent we can distinguish between various descriptions 
of the interior physical processes. This method can be applied to a wide range of 
astrophysical problems, and here I present one such example which shows that 
the convective core overshoot parameter can be constrained with one identified 
mode if the pulsating component is in an eclipsing binary system. 



1. Introduction 



It is of extreme benefit to have a reliable method to explore both the scientific 
potential and the limitations of a particular astrophysical system. This paper is 
dedicated to discussing how to use some properties of singular value decompo- 
sition (SVD) to do just this. We can use SVD to do such things as predict the 
precision of global quantities, explore the impact of extending a data set from 
3 months to 6 months (e.g. with Kepler), and probe the scientific benefit of 
complementary observables. If we understand the potential of the system, then 
we can carefully choose the scientific question that we would like to pursue. As 
a specific example, we will test if an eclipsing binary system containing a pulsat- 
ing component has sufficient information to allow us to constrain the convective 
core overshoot par ameter, a ov . This method of using S VD was exploited by 
Brown et al. ( 1994 ), and mo re re cently in the lite r ature, by Miglio &; Montalban 
(|2005l h lCreevev et al.1 (|2007l > and lCreevevI (|2008al h iPress et alJ (|l992h while con- 
centrating on SVD from a numerical point of view, also briefly discusses some 
of these properties. 



2. Observational data and interpretation 

Normally in astrophysical problems, one is confronted with comparing a set 
of observational data with an astrophysical model for scientific interpretation. 
This model should describe as best as possible the physical or chemical processes 
observed. The input to this model are a set of parameters of the system, which 
we shall denote by P. Often we know what the parameters are, but we do not 
know their quantities. For example, in a spectroscopic binary system, such 
parameters are the orbital period, the mass ratio and the inclination of the 
system with respect to the observer. Given a specific set of input P, this model 
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will produce a set of expected observables, which we denote by B. In the binary 
system example, B could be the radial velocities of both components. Scientific 
interpretation often then involves comparing the real observational data, O, to 
B to infer the values of P. 

If the system is linear, or if it can be described locally as linear, the follow- 
ing equation can determine the true parameters of the system, by calculating 
(iteratively) the parameter changes <5P to make to an initial guess Po: 

Q 13 

SP = VW _1 U — - — , (1) 

where e are the measurement errors. Here, we can see that the parameter 
changes needed to match the observations come from mainly a product of two 
terms: the latter term in Eq. [1] is simply the scaled difference between the 
observations and the expected observables — the discrepancies — and the first 
term in Eq. [1] is the inverse of the sensitivity matrix (calculated from stellar 
models), expressed in its SVD form. (This matrix is: ^p/e •) The amount by 
which we need to change any parameter to successfully match O with B comes 
from a product of the discrepancies in the observations, and a set of linear 
vectors (given by U and V) describing the relationship between each observable 
and each parameter. The rest of this paper is not dedicated to solving these 
equations, but rather to exploring the information contained in U and V to 
gain an understanding of the system under study. 

2.1. Relationship vectors versus the sensitivity matrix 

Let's concentrate on a more concrete example. Imagine that the system under 
study is a single isolated pulsating star. The model consists of a stellar evolution 
and structure model, coupled t o an adiabatic oscillation code. He re I use the 
Aarhus Stellar Evolution Code ( Christensen-Dalsgaardl[i~982l l2007al ) , coupled to 



the ADIabatic PuLSation code (|Christensen-Dalsgaardll2007bl ) . The P are mass 



M, age r, initial mass fractions of the elements X and Z, and a mixing-length 
parameter to describe the outer convection zone a. The B are the classical 
observables of effective temperature T e g, luminosity L* and metallicity [M/H], 
and the average seismic quantities: the small and large frequency separations, 
(Sis) and (Au). These quantities are the average values of the observed 8v^ n = 
u i,n ~ u i+2,n-i an d Af/ n = v\ n — v\ n —\ where v\ n are the oscillation modes of 
degree I and radial order n. The errors on each of the observations are typical: 
200 K, 0.1 L , 0.1, 1.3//Hz, 1.3/xHz. 

Table Q] shows an example of the sensitivity matrix of the system just de- 
scribed, taking as reference a star with parameters similar to the Sun. Each 
element is a partial derivative of one observable (specified in the leftmost col- 
umn) with respect to one parameter (specified in the top row). For example, 
the partial derivative of with respect to r is 50.5. Note that these values take 
e into account, and will change if we assume some other values of e. 

If there were a discrepancy in only one of the observables, say, (Sv), it is 
not so clear by reading this matrix which parameters and by how much each of 
these will need to change to reconcile O with B, mainly because (5u) is sensitive 
to each parameter. However, representing this matrix in its SVD form allows us 
to do this very simply. 
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Table 1. Sensitivity matrix: partial derivatives of some stellar observables 
with respect to the global parameters, divided by the measurement errors 



M r X Z a 

T cS 36.5 0.3 -70.5 -421.9 3.4 

58.5 50.5 -71.1 -390.0 6.2 

log(Z/X) 0.0 0.0 -6.1 217.3 0.0 

(Ai/) -3180.9 -58.1 3913.82 18874.3 235.5 

(5u) -264.0 -13.9 384.5 1628.7 4.18 
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Figure 1. The decomposition matrices U and V, describing the relationship 
between the stellar observables (left panel) and the stellar parameters (right 
panel) for a single 1 M Q star. 

Figure Q] is a graphical representation of the SVD of the matrix given in 
Table [TJ In the decomposition both U and V are orthonormal matrices that 
span, respectively, the observable and the parameter spaces. Therefore each 
element is a value that varies between -1 and 1, and in Figure [T] each value is 
represented by a triangle, whose magnitude and direction is proportional to it. 

Instead of thinking of these matrices as two separate matrices, we should 
interpret them rather as sets of vectors that are related according to their posi- 
tion in the matrix. For example, row 1 in U (vector 1 or U\) corresponds only 
to vector 1 in V and it is related by the first (highest) singular value W\. 

Now, we can investigate which of the parameters should be changed in 
order to reconcile the discrepancy in the observable (Sv). Inspecting the column 
containing this observable on the left panel, we see that the largest component is 
contained in the 5th row of the matrix U (t/5). This implies that the parameter 
changes given by V5 are those necessary to reconcile the discrepancy in the 
observation. However, V5 (the bottom row on the right panel) has mainly one 
non-zero value, and this corresponds to r. This means that in order to resolve 
the discrepancy in {8v), one would need to adjust the value of r. The amount by 
which r needs to be changed is proportional to the discrepancy in the observation 
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and inversely proportional to the corresponding singular value, in this case I/W5. 
W5 is the smallest singular value, so I/W5 is the largest inverse. Now, we can 
also see that resolving discrepancies in the observables that appear in vectors 
U±, U2, etc. would cause smaller SP than those in U5 because 1/W±, I/W2, etc. 
are much smaller than I/W5. This implies that the parameters that appear in 
the topmost vectors are only allowed to be adjusted by a small amount, while 
those in the lowest vectors can vary more, i.e. the parameters appearing in 
V\_, V2, etc. have tighter constraints, and thus have the smallest uncertainties 
associated with them. In fact, the uncertainties <jj for each j parameter come 
in a neat and compact form: 



Lets take as another example the observable (At/). This observables appears 
mainly in Ui, implying that the parameter adjustments given by V\ will resolve 
any discrepancy in this observation (if one exists). In this case, we would need 
to increase M and decrease X by the same amount, while also increasing the 
value of Z by a larger amount. Again, the actual value by which we need to 
adjust these parameters is proportional to \jW\ and the discrepancy in (Ais). 
We can now begin to understand that these matrices show quite directly which 
observables contribute to determining each parameter. 

As a third example, if the uncertainty in the mixing-length parameter a 
were quite large and it was in our interest to decrease it, given that a appears 
mainly in V3 and V4, a reduction in the errors on the observables appearing in U3 
and U4 would cause the desired decrease in a(a) — these responsible observables 
are mainly T e g and L* (in this particular case). 

One must also take some precaution with this interpretation: these results 
are sensitive to the observational errors that we assume, the set of observables 
that we take into account, and the range of parameters that we are studying. 
Here we assume a single star with solar characteristics, but these relationships 
will change if the star is at a different evolutionary stage, metallicity, and mass — 
for a given reference set of parameters, these relationships extend to parameter 
ranges where we can still assume linearity with respect to the reference set. 



I showed briefly how we can interpret the decomposition matrices, and how useful 
the properties of SVD are for understanding the system under study. Not only 
can we investigate the role that each observable plays in determining the system 
parameters, but also SVD shows us if an observable has no important role, 
i.e. if the observable is redundant (in the case of an over-determined system). 
Being redundant implies that this observable can be used to test specific physical 
phenomena because they are independent of the stellar model chosen to represent 
the star. 
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Figure 2. Parameter adjustments dP needed to most adequately fit O when 
the model used for the inversion is incorrect. The last four parameters are 
rotational velocity of component A, component B, the distance to the system 
and the inclination. The values for va and vb fall outside of the figure. 



3. Determining the convective core overshoot parameter 

Let's suppose now that the system is a main sequence eclipsing spectroscopic 
binary system with component masses of roughly 1.8 and 1.1 Mq and the 1.8 M Q 
is a 5 Scuti star. If we consider both photometric and spectroscopic light curve 
observations, then the observables will comprise of things such as the effective 
temperatures T e g_A, and ratio T e g_B/T e g_A, inclination i, orbital period IT, radii 
Ra, Rbi mass estimates M^sini, Mb sin i, and an identified oscillation mode 
v from the 5 Scuti component. An analysis of these observables following the 
discussion in Sec. 2 shows that all of the parameters of the s ystem (now incl uding 
two stars that we assume are coeval) are well-constrained ( Creevevll20 08b). We 



pose the following question: Can the single identified oscillation frequency be 
used to learn about the convective core overshoot parameter a ov ? 

In order to set about answering this question, we simulated a set of observa- 
tions O, from a model with a ov = 0.3. We then used EqfT] iteratively to recover 
the input set of parameters, while using the correct model in the inversion. The 
parameters converged to the original set to within a small amount, we denote 
these by Pf- We obtain the uncertainties, c(P) from Eq. [21 

Now we change the inversion model to a ov = 0.2, and we try to recover 
the original observations using again, Eq. [21 while use the fit Pf as the initial 
guess. Because we are using the incorrect physical model to fit the observations, 
then we should expect to see some non-zero values of SP. The question is, is the 
difference in the observables sufficient to be able to detect? We can argue that 
if there are <5P that are larger than the <x(P) calculated, then we conclude that 
there is enough information in B to detect an error in the physical assumptions. 

Figure [2] shows <5P/<x(P) for each P after the inversion with the incorrect 
(a = 0.2) model. The set of observables include one identified mode. The 
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Figure 3. Discrepancies in O when they are fit to the a ov = 0.2 model. 
The best set of P result in large deviations for some observables. 



horizontal dotted lines represent ±l-<x(P). It is clear that the eclipsing binary 
observables are capable of providing different Pp solutions based on diverging 
physical assumptions — implying that there is sufficient information there to 
learn about a ov . 

But how can we distinguish between the two fit parameter sets given that 
both make physical sense? Figure shows the new discrepancies in the observ- 
ables, where B are calculated from the a ov = 0.2 models (our assumption for 
the inversion) with the fit parameters corresponding to those in Fig. [2j The 
large discrepancy in three of the observables indicates that our model is incor- 
rect. The identified mode is the observable that provides most evidence that 
a ov = 0.2 is incorrect, its discrepancy is 40 times the allowed amount away from 
the observed value. Repeating this for various values of a ov allows us to recover 
correctly the input value of 0.3. 
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